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ABSTRACT: We use the recently conjectured exact S'-matrix of the massive O(n) 
model to derive its form factors and ground state energy. This information is then used 
in the limit n —>■ to obtain quantitative results for various universal properties of self- 
avoiding chains and loops. In particular, we give the first theoretical prediction of the 
amplitude ratio C/D which relates the mean square end-to-end distance of chains to the 
mean square radius of gyration of closed loops. This agrees with the results from lattice 
enumeration studies to within their errors, and gives strong support for the various as- 
sumptions which enter into the field theoretic derivation. In addition, we obtain results for 
the scaling function of the structure factor of long loops, and for various amplitude ratios 
measuring the shape of self-avoiding chains. These quantities are all related to moments 
of correlation functions which are evaluated as a sum over m-particle intermediate states 
in the corresponding field theory. We show that in almost all cases, the restriction to 
m < 2 gives results which are accurate to at least one part in 10^. This remarkable fact 
is traced to a softening of the m > 2 branch cuts relative to their behaviour based on 
phase space arguments alone, a result which follows from the threshold behaviour of the 
two-body S'-matrix, S(0) = —1. Since this is a general property of interacting 2d field 
theories, it suggests that similar approximations may well hold for other models. However, 
we also study the moments of the area of self-avoiding loops, and show that, in this case, 
the 2-particle approximation is not valid. 
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I. INTRODUCTION 



In the past few years, there has been remarkable progress in obtaining exact results 
for two-dimensional field theories 1 ' 2 . One application of these is to the scaling limit of 
classical statistical mechanics systems, close to a second order phase transition where the 
correlation length £ is much larger than the microscopic cut-off. Yet, besides checking 
the values of critical exponents and other universal finite-size scaling amplitudes which 
arise from studies of massless, or conformal, field theories, there has been little in the way 
of direct confrontation with statistical mechanics^. This is partly because the simplest 
predictions of massive quantum field theories are in terms of their particle content and 
S'-matrix elements, as quantum field theories in 1 + 1 dimensions. In order to make 
contact with easily observable quantities relevant to statistical mechanics, it is necessary 
to go further and derive off-shell behaviour. More recently, these problems have become 
tractable, at least for a class of integrable theories, with the development of methods 
for calculating off-shell form factors 4 ' 5 and the technique of the thermodynamic Bethe 
ansatz 6 ' 7 . 

However, another obstacle is the lack of precise information with which to compare 
such predictions. Aside from the Ising model 8 ' 9 little is known about off-critical corre- 
lation functions in statistical models. Ideally, of course, one should compare with real 
experiments, but unfortunately it is notoriously difficult to obtain reliable data on two- 
dimensional systems, uncontaminated by the effects of impurities and finite-size. While, 
in principle, it is possible to carry out high-resolution scattering experiments with X-rays, 
such studies appear to have gone out of fashion with experimentalists. 

We are therefore left with comparison with the results of numerical studies. In contrast 
with real experiments, these are actually easier to carry out for large systems in low dimen- 
sions. While Monte Carlo studies can study quite large systems, their results are subject to 
statistical errors which make comparison with detailed predictions of the correlation func- 
tions difficult. Since we wish to study such functions unaffected by finite-size effects, exact 
diagonalisation of finite transfer matrices is also not useful for our purpose. This leaves 
series expansions. The problem for which the longest such expansions have been carried 
out, and therefore the most accurate information is available, is that of self-avoiding walks 
(see Ref. 10 for a review of the current state of knowledge, and for further references). 
Although this is strictly not a problem in critical behaviour, its relation to the critical be- 
haviour of the O(n) model in the limit n — > was shown long ago by de Gennes 11 , and it is 
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possible to express all the quantities relevant to the asymptotic behaviour of iV-step walks 
in the limit N — > oo in terms of correlation functions and amplitudes. The correspondence 
actually relates the fixed fugacity ensemble, where each step is weighted by a factor x, 
to the n — > limit of the O(n) model, at a temperature simply related to x. One result 
of this is that the actual correlation length of the O(n) model is not directly measurable 
through self-avoiding walks. 

The field theory of the two-dimensional O(n) model has been studied extensively. While 
for values of n > 2 this model is always in its disordered phase at non-zero temperature, 
its analytic continuation to —2 < n < 2 may be carried out, and it does then have a finite 
T c . In fact, the lattice version of this model at criticality may be mapped via standard 
techniques onto the Coulomb gas model 12 ' 1 *', with the result that the critical exponents 
previously conjectured 14 ' 15 and the value of the central charge c in the corresponding 
conformal field theory may be easily read off. The standard methods of conformal field 
theory 1 ^' 1 ^ then show that the energy operator corresponds to a degenerate representation 
of the Virasoro algebra, labelled by its position (1,3) in the Kac table. It was shown by 
A. Zamolodchikov 1 ^ that a conformal field theory perturbed by a relevant operator of this 
type is integrable, in the sense that an infinite number of the conserved charges in the 
conformal field theory survive the perturbation. If the perturbed theory is massive, which 
we expect to be the case for T > T c in the O(n) model, the existence of these conserved 
charges implies that scattering in the theory is purely elastic, and that the many-body S- 
matrix factorises into a product of two-body scatterings 19 . For the O(n) model one expects 
to find a simple particle spectrum: an n-plet of particles states transforming according to 
the vector representation of O(n), and, since the interaction is repulsive, no bound states. 

On this basis, A. Zamolodchikov 20 conjectured an exact S-matrix for this theory, based 
on the principles of analyticity, unitarity, crossing, and the Yang-Baxter equations. His 
result represents the simplest solution to these conditions, but it is not unique, because of 
the CDD ambiguity. 

Once the S-matrix is assumed, however, there is a well-defined program for investigat- 
ing the off-shell behaviour of the theory, through the form factors 4 ' 5 . For a local field <f)(r), 
these are defined as matrix elements (O|0(O)|/?i, a\; . . . ; f3 m , a m ), where the ket represents 
an asymptotic m-particle state labelled by the rapidities (3j and O(n) colour indices cij of 
the particles. The form factors satisfy a set of equations relating them to the S'-matrix, and 
other requirements of analyticity and crossing. Once again, it is generally assumed that 
the simplest, or 'minimal', solution of these equations should be chosen. This program 
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has been carried through in detail only for theories with non-degenerate single-particle 
states 21 ' 22 ' 23 ' 24 , or with only simple symmetries 4 ' 5 ' 25 ' 26 . The spectrum of the 0(n) model 
is, however, more complicated, and the two-body .S-matrix acts on the direct product of 
two n-dimensional representations. While, for specific integer values of n, it is clearly 
possible to decompose this into irreducible amplitudes for which the S'-matrix is diago- 
nal, this is not useful for continuation to n = 0. We are therefore forced into additional 
complications in order to include the O(n) structure. 

Fortunately, in most cases this may be avoided for all intents and purposes. The 
purpose of computing the form factors is to reconstitute the two-point correlation functions 
through the unitarity sum 

<0(r)0(O)> = £ / \(0\cf ) \p 1 ,ar,...;P m ,a m )\ 2 e- Mr ^ coshpj . (1) 

m=0 colours ^P^ 80 s P ace 

Clearly, for the infrared behaviour (large r) only the lowest values of m are important, 
while, in the ultraviolet region of small r, in principle all the intermediate states contribute. 
However, for our purposes, we shall be interested in moments of the correlation functions 
of the form J r p {4>{r)4>(Q))d?T. For large p, the large r end of the integration region will 
be emphasised, and the truncation to, say, m < 2 should be valid. Remarkably, however, 
we find that, this can be true for p = 2 and even p = 0. In the case when is the energy 
operator (which couples only to states with m even), the value of the p = 2 moment is 
known exactly from a sum rule derived from the c-theorem 27 ' 28 . The result obtained from 
the truncation to two-particle states agrees with this up to about one part in 10 3 . This 
approximation is less good for the p = moment (having an error of about 10%), but, as 
will be discussed below, this moment may be found exactly from the TBA. When is the 
O(n) vector field itself (corresponding to the spin of the lattice O(n) model), it couples 
only to states with odd m. Then the truncation to m = 1 appears to be very accurate 
even for the zeroth order moment. 

The reason for this remarkable numerical simplification appears to be that the form 
factors to the states with higher numbers of particles have a much softer behaviour close 
to the m-particle threshold than expected on the basis of phase space alone. In fact, 
each form factor contains an explicit factor of Yli<j(Pi ~ Pj)- This may be traced to the 
threshold behaviour of the two-body S- matrix, S(0) = — 1. In general, unitarity implies 
that S(f3)S(—f3) = 1, so that, in principle, S(0) could be ±1. However, the upper sign 
appears to be realised only in a free boson theory: otherwise, in all interacting theories, 
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the lower sign holds. Thus, this suppression of the higher particle states would appear to 
be a very general feature of two-dimensional theories. 

In the case of the O(n) model, the restriction to m < 2 makes the otherwise cum- 
bersome problem of the group theory factors straightforward. As a result, we are able 
to obtain accurate values for the moments of the spin-spin and energy-energy correlation 
functions in the n — > limit, based, of course, on all the assumptions of minimality for 
the S-matrix and the form factors referred to above. As will be shown in detail in the 
subsequent sections, these moments are simply related, respectively, to moments of the end- 
to-end distance of self-avoiding walks, and of the mean square distance between points of 
self-avoiding loops. For these lattice problems, extensive numerical work has been carried 
out on the mean square end-to-end distance of A-step walks, (i?e)iV ~ CN^ U ^, and the 
mean square radius of gyration of loops, (R%)n ~ DN^ V ^®, where, from Coulomb gas ar- 
guments, v = | 12 , and C and D are amplitudes. By themselves, they are not universal^ 
(as will also follow from our analysis below), but the ratio C/D should be free of all metric 
factors and therefore universal. 

The only quantity entering the calculation of this ratio which is not determined to high 
accuracy by the 2-particle truncation of the form factor approach is the zeroth moment Uq 
of the energy-energy correlations. But this is just proportional to the specific heat, and 
therefore may be found in an alternative way by differentiating the extensive part of the 
free energy. This may be derived from the thermodynamic Bethe ansatz. The TBA has 
been applied to the unitary minimal models with c = 1 — 6/(fc + l)(k + 2), with k integer, 
perturbed by the (1,3) operator, by Al. Zamolodchikov 7 . Although the O(n) model is 
not identical in its operator content to the minimal models, the results for the free energy, 
when perturbed by the same type of operator, must be the same. Thus we may take over 
Zamolodchikov's result, the only generalisation necessary being the continuation to non- 
integer fc, so as to allow the n — » limit to be taken. This turns out to be straightforward, 
and leads to an exact result for the universal part of the free energy per correlation volume, 
which is related to Uq. 

Piecing together all this information, we find C/D ~ 13.70. This is to be compared 
with the estimates 1 ^ 1 of 13.69 for the square lattice, and 13.72 for the triangular lattice. 
Thus the errors from the 2-particle truncation are less than the systematic errors of the 
current data from enumerations. This agreement is a non-trivial test of the O(n) model 
S'-matrix and form factors. (By contrast, the ratio C/D for ordinary random walks is 12.) 

The mean square radius of gyration of loops is related to the ratio of the p = 2 and 
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p = moments of the energy-energy correlation function. The generating function for 
all the moments is proportional to the mean structure factor of iV-step loops, S(q) = 
^2iji eiq '( ri ~ rj ^), which could be measured, for example, by light scattering. Since 
the 2-particle approximation gives the zeroth and second moments to such remarkable 
accuracy, and its accuracy is supposed to improve for the higher moments, this gives a 
very precise way of computing S(q), at least for moderate values of its scaling argument 
q 2 (R 2 g ). Thus we are able to make the first accurate calculation of this function. The higher 
moments are dominated by the behaviour of the 2-particle cut near threshold. Since, as 
we remarked above, this is softened due to the fact that S(0) = — 1, these moments behave 
as a function of p in a very different way from those of a free theory. Such a free theory 
describes random self-intersecting loops. Thus we are able to picture in a rather direct 
way how the condition of self-avoidance influences the statistics of the loops. On the 
other hand, the behaviour of S(q) at large values of its scaling variable qR is proportional 
to (qR)~^l v , as expected for an object of fractal dimension 1/v. The coefficient of this 
behaviour may be determined by comparing conformal perturbation theory with the results 
of the thermodynamic Bethe ansatz. 

Although the truncation to states with m < 2 works remarkably well for moments of 
the spin-spin and energy-energy correlations, this is not true for all correlation functions. 
For example, the O(n) model possesses conserved currents whose integrals generate the 
O(n) symmetry. It turns out that moments of the current-current correlation function 
are related to moments of the area (a p }]\[ of iV-step self-avoiding loops, in the limit n — > 
10 ' 31 . These are expected to scale as Thus we might expect to be able to 

compute the amplitudes in the 2-particle approximation. The amplitudes of both the 
first and second moments have been estimated from enumerations. However, it turns out 
that for these values of p, the 2-particle approximation is hopeless. This may be traced 
to two causes: first, current conservation forces an additional softening of the 2-particle 
contribution to the current-current correlation function near threshold; and second, the 
fact that the current has unit scaling dimension means that the ultraviolet behaviour is 
more singular, and gives a large contribution to the lower moments which is not adequately 
represented by the 2-particle contribution. 

While the form factor approach has, so far, been applied to the computation of two- 
point functions, it is not restricted to this case. Several amplitudes governing the statistics 
of open self-avoiding walks, for example their mean square radius of gyration, are related to 
moments of higher-point correlation functions. We show how in principle this calculation 
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may be carried out, and that, in the simplest approximation, we get results which agree 
with those of numerical lattice calculations to within a few percent. 

The lay-out of this paper is as follows. In the next section, we make precise the 
connection between the self-avoiding walk problem and the O(n) model, and in particular 
relate the various amplitudes and scaling functions whose values we are trying to compute 

to moments of appropriate correlation functions in the field theory. In Sec. (3) we recall 

on 

A. Zamolodchikov's arguments for the S'-matrix of this theory, and derive the 2-particle 
form factors. We show how this leads to an accurate estimate for the second moment of 
the energy-energy correlation function. In the next section we derive the zeroth moment 
of this correlation function by way of the thermodynamic Bethe ansatz, generalizing the 
arguments of Al. Zamolodchikov^ 7 so that the limit n — * may be taken. This result, 
combined with those of the previous section, then leads to our estimate for C/D. 

In Sec. (5) we present the calculation of the structure factor and comment on its fea- 
tures, and in Sec. (6) we discuss the current- current correlation function and show why the 
2-particle approximation fails for its lower moments. In Sec. (7) we consider the radius of 
gyration of open chains, which involves the calculation of a three-point correlation func- 
tion, and, finally, we summarise our conclusions and discuss possible further extensions of 
our results. 



In this section we summarize the well-known relation between these two problems in 
the limit n — > 0, in order to make precise the quantities we shall need later. 

Consider first the set of A-step self-avoiding walks on a given lattice. It is convenient 
to consider lattices of co-ordination number three (e.g. a honeycomb lattice in two di- 
mensions), but, later on, an appeal to universality will imply that in the critical region 
the results are universal. In the same spirit let us consider a specific lattice model with 
O(n) symmetry, in which the degrees of freedom are spins s a (r) at each site r, labelled 
by a 'colour' index a which takes values from 1 to n. The spins are normalised so that 
Tr s a S5 = 5 a 5, and the partition function is 



2. THE 0(A) MODEL AND SELF-AVOIDING WALKS 
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where the product is over all nearest neighbour pairs of sites (rj, rj). Once again, we expect, 
on the grounds of universality, that the particular form of the hamiltonian 
Hi = — Xlnn m (l + x s a( r i) s a(^j)) implied by (2) should not be important in the criti- 
cal region. In particular, we could consider the more standard hamiltonian 

H 2 = -z£ nn Ea s a( r *) s a( r j)- 

The continuum limit of either of these models close their respective critical points 
is assumed to be described by an 0(n)-invariant field theory, and, at the critical point, 
by a conformal field theory whose central charge c(n) = 1 — 6/(fc + l)(k + 2), where 
n = 2 cos(7r/(/c + 1)) with k > 16 ' 17 . We note, in particular, that c vanishes when n = 0, 
as it should, since it measures the finite-size correction to the free energy which would 
vanish in this case. Its derivative dc/dn at n = will be important: this takes the value 
5/37T. The scaling dimension of the energy operator is x e = 2k/ (k + 2) and equals 2/3 
for n = 0. It is related to the exponent v characterising the divergence of the correlation 
length by v = 1/(2 — x e ) = 3/4. This also fixes the critical index a = 2 — 2u of the 
singular part of the free energy. The other exponents which will enter are 7 = 43/32 which 
gives the divergent behaviour (x c — x)~^ of the susceptibility in the limit n = 0, and the 
magnetic scaling index rj = 5/24. They are related by the scaling equation 7 = z/(2 — r/). 

On expanding the product in (2), it may be written as a sum over self-avoiding loops: 

g _ ^ ^ ^.number of links ^number of loops ^ 

loop 
configurations 

As n — > 0, only configurations with a single loop survive in the 0(n) term. Thus, if is 
defined as the number of iV-step self-avoiding loops per lattice site, then 

NsJ^Pnx N = lim n _1 lnZ , (4) 

where J\f s is the total number of lattice sites. Note that the right hand side is proportional 
to the free energy, so is also proportional to the total area A. The free energy of the O(n) 
model is believed to be an analytic function of its temperature-like variable x for sufficiently 
small x, and to have a singularity of the form (x c — x)^~ a at some finite, lattice-dependent 
x c . Since the coefficients of the series on the left hand side are all non-negative, the 
singularities of the free energy on the circle \x\ = x c will determine the behaviour of pjy at 
large N . Hyperscaling implies that, in two dimensions, A~^~ InZ ~ nU£~^, where £ is the 
correlation length and U (which we shall calculate exactly in Sec. (4)) is universal. The 
correlation length itself has the critical behaviour 

t~Z (l- x / Xc )-« , (5) 
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where v = 3/4 but the metric factor £q is no n- universal. Putting all these factors together 
we see that 

J2pn* N ~aot 2 U(l-x/x c ) 2v , (6) 
N 

as x — > x c , where ag = A/Af s is the area per site. This implies that, as N — > oo, 
~ BN~ 2l> ~ 1 fj, N , where = a^T 1 and the the amplitude £? is given by 

For certain lattices p^y actually vanishes by symmetry except when iV is divisible by an 
integer a. This implies that the left hand side of (6) is actually a series in x a : and therefore 
has a equivalent singularities on its circle of convergence, thus accounting for the lattice- 
dependent factor of a on the right hand side of (7). For example, for the square lattice 
a = 2. 

Next, we consider the energy-energy correlation function. Close to the critical point, 
we write H 2 = H% + (x c - x) Ebonds r £ lat ( r )> where £ lat (r) = ^ a s a(r<)s a (r > ) is the 
lattice energy located on the link r which connects the sites r< and r>. The correlation 
function (£^ at (ri) £^' At {r 2 )) then receives contributions from all self-avoiding loops which 
contain the links r\ and r 2 . If the number of such loops with N steps is /_/v( r l? r 2) 5 then 

nJ2fN(n,r 2 )x N - 2 = (£ lat (ri) £ kt (r 2 )> • (8) 
N 

In particular, if we sum over all v\ ^ r 2 on the left hand side, we sum over all such iV-step 
loops, and obtain 

nAf s J2N(N-l)p N x N - 2 = J^^in)^^)} , (9) 

which is just the second derivative of (4) with respect to x. The right hand side is just the 
specific heat of the O(n) model. 

In the same way, if we multiply by (r\ — r 2 ) 2 before summing, we get the generating 
function for weighted by 

2 

J](ri-r 2 ) 2 = (Vl - r) - (r 2 - f)) =2N 2 R 2 g , (10) 

ri,r 2 ri,r 2 

so that 

2uKsJ2 n2 Pn(RI)nx N - 2 = J2( r l-r 2 ) 2 (£ l&t (n) £ lat (r 2 )> . (11) 

N n,r 2 
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Similarly, the generating function for the structure factor Sjy(q) = N 1 J2 ri r 2 ( e * 9 ^ T2 ^) 
is given by 

n^ s J2N PN S N (q)x N - 2 = e^- r ^(S^(n) £ lat (r 2 )> . (12) 
N n,r 2 

In the scaling region, the dominant contribution to the sum on the right hand side comes 
from \r\ — r 2 \ of the order of the correlation length, and we may therefore use a continuum 
field theory description. The lattice energy £ lat is replaced by the continuum energy density 
£{r) in such a way that they enter into the hamiltonian (action) in the same way: 

J~V at (r) -> j 8{r)Sr . (13) 
r •* 

However, in the field theory, the energy operator is special in that it represents the per- 
turbation (x c — x) J £{r) (fir of the critical theory. It is therefore proportional to the trace 
0(r) = Tjf of the stress tensor. Explicitly 

0(r) = -2nv- l (x c -x)E(r) , (14) 

where the conventional factor of 2tt has been included in the definition of the stress tensor. 
Since the normalisation of the stress tensor is fixed, its correlation functions are completely 
universal. In particular 

<e(r)e(0))=<- 4 <h(r/6 , (15) 
where $i is a universal scaling function. 

From (11) we then see that, in terms of this scaling function, 

2n n2 Pn(R 2 ) N% N ~ 2 = ao(u/2n) 2 (x c - x)' 2 j r 2 {®{r)Q{Q))d 2 r (16) 

= nao(v/2Tr) 2 U 2 (x c -x)- 2 (17) 

where U 2 = j p 2 <&i(p)cfip is universal, so that 

2Np N {R 2 g ) N ~< T a (v/27r) 2 U 2f x N . (18) 

The right hand side of (16) is in fact given by a sum rule which is a consequence of 
Zamolodchikov's c-theorem 27 : 

c(n) = (3/4tt) J r 2 (O(r)O(0))d 2 r , (19) 
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where c(n) is the central charge of the 0(n) model. In terms of the amplitude D in the 



relation (Ft*) ~ DN 2v , this implies the relation 



SD= 3272 aa ° ' (20) 

which was derived in and generalised in Alternatively we see that the c-theorem 
sum rule gives an exact value for U<i'- 

U 2 = J p 2 $l(p)d 2 p=™ • (21) 

From (9), the amplitude B may itself be written in terms of the zeroth moment of $1, 
since 

J2 N 2 PN x N - 2 = a (is/2n) 2 (x c - x)~ 2 j (0(r)0(O))d 2 r = a (v/2n) 2 (x c - x)-\- 2 U Q 

where Uq = f <&\(p)d 2 p. Hence we obtain the alternate result for the amplitude B: 

B = ™0t0 2 (»/^) 2 T{2 U ° 2u) ■ ( 22 ) 
Comparing with (7), this gives the further sum rule 

U = J <$> l {p)d 2 p = (2tv/u) 2 (-2u)(1-2u)U . (23) 

We stress that this is merely a consequence of the fluctuation sum in (9) being proportional 
to the specific heat. As we shall see in Sec. (4), the amplitude U may be obtained exactly. 
The above two sum rules for Uq and U<i then give an estimate of the errors of the two- 
particle approximation to (0(r)0(O)). 

The amplitude D defined above for the radius of gyration may also be written, from 
(7, 18, 23) in the more intuitive form 

n _ c2 r H U 2 (24) 

D - Co 2i> + 2v) TT ' (24) 

As expected, the mean square radius of gyration is proportional to the ratio of the second 
to the zeroth moments of the energy-energy correlation function. 

Finally, the structure factor is related to the Fourier transform of $i by 

E Np N S N ( q )x N ~ 2 = a (u/2n) 2 (x c - x) 2 C 2 [ e^* $i(p)d 2 p . (25) 
N J 
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Disentangling the large N dependence from this equation is more complicated, and will 
discussed in Sec. (5). 

Next, we consider the case of linear self-avoiding walks. Let cjy(r) be the number of 
such walks with N steps from the origin to the site r. The generating function for this is 
just the spin-spin correlation function of the O(n) model as n — > 0: 

J2 c N(r)x N = lim( Sl (r) Sl (0)> . (26) 

In the scaling region, we expect this correlation function to scale as 

(si(r)ai(O)) ~ ftT^Cr/O > ( 27 ) 

where $2 is universal, but the metric factor b is not. Note that the magnetisation is 
different from the energy operator in that it has no absolute normalisation. From (26) we 
may easily compute the generating function for the number c/v of all iV-step walks: 

c N x N = b^-\l - x/x c )-^ V , (28) 

N 

and for this number weighted by their squared end-to-end distance R\: 

J2 c N(R 2 e)NX N = b$-\l - x/x c )-^V 2 , (29) 
N 

where Vq and V2 are the zeroth and second moments, respectively, of $2- We then find, 
for the amplitude C in the asymptotic behaviour ~ CN 2iy , the result 

C -^fWT^)V ■ (30) 

Note the close similarity with (24). (The extra factor of 2 in the denominator in (24) is 
due to the factor of 2 appearing in (10).) 

Neither of the amplitudes C or D is universal, because of the metric factor £q. However, 
this cancels in their ratio, and we find 

1 r( 7 + 2^)r(«) v u 2 ■ {6L) 
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3. ^-MATRIX AND FORM FACTORS OF THE 
THERMAL PERTURBATION OF THE 0(N) MODEL 



3.1 Scattering theory 

Since in the scaling limit the energy operator of the O(n) model on the lattice corre- 
sponds to the primary operator </?i 3 of the conformal model, the thermal perturbation of 
the critical point action is described by 

A = A c + (x c -x) J ip 1:3 (x)d 2 x . (32) 

As with any ^3 deformation of the minimal conformal models, the QFT defined by the ac- 
tion (32) presents a number of higher integrals of motion which guarantee its integrability 18 . 
For x < x c (the only case that we will investigate) the model develops a finite correlation 
length and the associate scaling QFT is mainly characterised by the factorised S-matrix 
of its massive excitations. The scattering theory has been proposed by Zamolodchikov 20 
and its main features may be summarised as follows 1 . 

First of all, on the basis of the form of the partition function (4), Zamolodchikov argued 
that it is possible to interpret the loops as trajectories of a set of n particles Ai{(3) (with 
mass M) that belong to the vector representation of O(n). Hence the scattering matrix for 
the process |^4j 1 (/?i)^4j 2 (/?2)) — * |^4j 1 (/5i)^4j 2 (/?2)) may be written as 

SilSiM = So(Pu) + ^(A2) t 5 h + ^12) h^ in , (33) 

where flyx = f3\ —fa- So is the amplitude for the transmission process whereas <Si and S2 are 
respectively the amplitudes for the reflection and annihilation processes. Equivalently, we 
may decompose the scattering matrix into channels of definite isospin: for the symmetric 
traceless, antisymmetric and isosinglet channels we have respectively 

S S ((3) = S ((3) + sm 

S A ((3) = S Q (P) - Sm (34) 
Sm = S Q {(3) + SM + nS 2 (P) • 

It is worth mentioning that the same S'-matrix, although with a different interpretation, 
has been also derived by Smirnov 33 using a quantum group reduction of the scattering 
theory of the Sine-Gordon model. 
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The amplitudes S\ and £2 are linked to each other by the crossing symmetry relation 



Sl(P) = 5 2 (t7T - P) , 



(35) 



whereas Sq is a crossing symmetric function. 

In order to take into account the property that the loops entering (4) are non-intersecting 
paths, Zamolodchikov suggested imposing the condition 



S O {0) = 



(36) 



Using the Yang-Baxter equations and the unitarity condition, the final form of the minimal 
.S-matrix is then given by 



Si(P) = — sinh 
S 2 {j3) = - sinh 



in 







k+1 

g 
k+1 



R{(3) 



(37) 



where 



R(P) = 



1 



exp 



DO 



dx 



sinh sin xft 



x sinh iih+lh cosh if 



Its analytic structure may be read off from its infinite product representation 



(38) 



R(0) 



1 



( 1 ft. "\ 00 F 

\ L m(k+l)) -i-r 1 



( Jl §_ 



Sm7T 



1 



ft 



r 1 + 



k+1 in(k+l) 
21 ft 



r 1 + 



X 



k+1 m(k+l) 



1 = 1 1 

k+1 ^ iw(k+l) 



\k+l m(k+ 



it) 



in(k+l) 



21 



+ 



ft 



r[+^- + 



k+1 ^ in(k+l) 

I [1 + ++T + 



ft 



k+1 ^ iTT(k+l) 



(39) 



r 1 + 



21 



+ 71 



k+l^ in(k+l)J L \k+l in (k+1) J ^ ^ ^ k+1 in (k+1) t 

Since there are no poles in the physical sheet < Im/3 < tt, no bound states appear. 
Hence the whole particle content of the model only consists in the n degenerate states 
of the vector representation of O(n). For n = 1 the S'-matrix of the isosinglet channel 
correctly coincides with the S'-matrix of the thermal perturbation of the Ising model 



ft 



r 1 + 



21-1 



ft 



whereas in the limit n 



SiW) + S 2 {(3) = -1 , 

0, the final form of the scattering amplitudes is given by 

S 1 (P)=-G(P) 
S 2 (P)=-itanh(^\G(P) , 



(40) 
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where 



G((3) = exp 



00 dx sinh^sinx/? 
x sinh Tlx cosh ^ 



(41) 



Notice that the S'-matrix of the 0(n) model possesses nontrivial asymptotic behaviour for 

P -> ±00 



,±zttA 



where 



3A; + 2 



(42) 



2(fc + 1) 

According to 26 , this implies that the particle operators of the O(n) model satisfy a gener- 
alised statistics. 



3.2 Form Factors 



The previous discussion of the scattering theory of the thermal deformation of the O(n) 
model gives rise to a series of questions related to the minimality assumptions made in its 
derivation or to the non-integer values assumed by the variable n. In addition the final 
scattering amplitudes are not directly related to the underlying microscopic formulation 
of the model, making it apparently problematic to judge their validity. To answer such 
questions, it is useful to recall that for integrable models the knowledge of the S- matrix is 
a powerful starting point to compute their correlation functions through the form factor 
bootstrap approach 4 ' 5 . Hence we may take the point of view of considering the O(n) 
scattering theory discussed above as simply the basic tool to implement the program of 
reconstruction of correlation function for this model, the whole justification of the approach 
being in the final comparison with quantities extracted by other methods. 

In the form factor bootstrap approach the correlation functions are computed by ex- 
ploiting their spectral representations, that is, their expression as an infinite series over 
multi-particle intermediate states. For instance, the two-point function of an isoscalar 
operator <&(x) in real Euclidean space is given by 

OO 

(*(s)*(0)) = V f dPl U0 • Z n < Ol$(a:)lft, 01; ..p n , a n >< Pi,a v .p n , On|$(0)|0 > 
^ Q J n!(27r)" 

= E/ f? |^i<...,«„(A---/'»)! 2 -p(-^E<« h ft) . 

n=0 V i=l I 

(43) 
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where a sum on the colour indices is implied. In (43) r denotes the radial distance, i.e. 
r = \J Xq + x\ and 

. . . , P n ) ai ,a 2 ,...,a n = < | $(0) | A, a i; . . . ; n , a n > (44) 

are the so-called form factors (FF). The normalization of the asymptotic states is fixed as 

(ft, ai I 02, a 2 ) = 27r<y ai , a2 <y(0i -0 2 ) • (45) 

Since the spectral representations are based only on the completeness of the asymptotic 
states, they are general expressions for any QFT. However, for integrable models, they 
become quite effective because the exact computation of the form factors reduces to finding 

a solution of a finite set of functional equations. In fact they satisfy the Watson equations 

4,5 

F (Pi...Pi+i,Pi, ■■■Pn) ai ..at +1 ,a' i ..a n =F (01"0i> Pi+l~Pn)a 1 ..a i ,a i+1 ..a n ^a'^C+S^ ~~ 
F®([3i, . . . , p n + 2m) ai ,...,a n =F®(P n , • • • , Pn-\)a n ,a u -,a n -i • 

(46) 

The first of (46) show that the monodromy properties of the FF are determined by the 
two-body S'-matrix of the model. On the other hand the second equation states that 
an analytic continuation in the variables Pi simply induces a reordering in the set of the 
asymptotic particles. Notice that this system of equations does not uniquely fixed the 
solution since the multiplication of F®(Pi, . . . , Pn)ai,...,a n by a symmetric, 2iri periodic 
function leaves (46) untouched. For a specific operator this ambiguity may be generally 
solved by the knowledge of the asymptotic behaviour of its form factors and their analytic 
structure. Let us briefly discuss the two aspects separately. 

The asymptotic behaviour of a FF under a simultaneous shift of all rapidity variables 
is simply dictated by relativistic invariance 

F $ (/3i + A, 02 + A, • • • , Pn + A) ai) ...,a„ = e sA F $ (/3i, ft, ... , 0n)a 1 ,...,a„ , (47) 

where s is the spin of the operator Notice that for scalar operators the FF depend 
only on the differences 0^ = 0^ — Pj. Secondly, restricting our attention only to those 
operators which have two-point functions with a power-law ultraviolet behaviour, we have 
to require that their FF behave asymptotically no worse than exp(/c/%) for 0^ — > oo, where 
k is a constant independent of i. If a perturbative formulation of the theory is available, 
the constant k may be fixed by matching the asymptotic behaviour of the FF with the 
asymptotic behaviour of the relevant Feynman diagrams which contribute to the amplitude. 
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Concerning the analytic nature of the form factors, their pole structure gives rise to 
a set of recursive equations which relate the /-particle FF (/ > 2). Since in the O(n) 
model there are no bound states, the only singularities which appear in the FF are those 
associated with the annihilation poles [3 a = [3^ + in with residues 5 



2n i resF $ (/3 + in, (3, fa, . . . , Pi)a ,a ,a u ...,ai = 

(« 2 ■ ■ ■ si - Sai;:::i(Pi, • • • , a i /?)) . . . , a)^,...^ 



where 



(48) 



sfci^i, ...,Pi\p) = <\'£(/?i - 0)S%£{/32 -P)... s a a l;ZM - P) ■ (49) 



These equations induce a recursive structure in the space of the FF relating matrix elements 
with / + 2 and / particles. 

After this short discussion on the general properties of the FF, let us turn our attention 
to the O(n) models. In order to analyse the correlation functions in the high-temperature 
phase of the model, it becomes convenient to define initially a two-particle form factor 
Fmin(P) that does not depend on colour indices, does not contain poles on the physical 
sheet and has the mildest behaviour for large values of f3. The first requirement implies 
that the monodromy property of such FF is induced by the S-matrix of the isosinglet 
channel Sj(f3) i.e. 



F m M = [Sm + nS 2 (P)]F min (-P) 
F min (iir-P) = F min (i7r + P) . 
The solution of these equations which satisfies the above requirements is given by 



(50) 



,(13) = sinh^ exp 



dx 



■ u irkx ■ 2 fix 
sinn !L y ± sin tj- 



x sinh Tix sinh 7r ( fc + 1 ) x cogh 



7TX 

T 



(51) 



where (3 = in — (3 and 



Skip) 



sin(/3/2) 



(52) 



sin (P/(k + 1)) 

Note that f\{f3) = 1. In the self-avoiding case (n — * 0) the monodromy property of F m [ n (/3) 
is dictated only by Si((3) of (40) and its expression reduces to 



mm 



((3) 



■ ^ 

sinn — exp 



'" 50 dx sinh ^ sin 2 \(in - (3)x 



x 



sinh 2 nx cosh ^ 



(53) 



Useful properties of this function are discussed in the Appendix. 
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Notice that at the threshold ((3 — > 0) F m [ n ((3) goes linearly to zero since 5/(0) = — 1. 
The vanishing of the F m [ n (f3) at the threshold has far-reaching consequences because it 
induces a suppression of the higher particle contributions in the spectral representation of 
the correlation functions. In order to see this, let us consider the general parametrisation 
of a FF of an operator $ a v,"z w ith colour indices a, . . . , m 



< | ^'-'-(0) | ft, ai ; . . . ; Pi, ai >= tf-Xifh, ,Pl)U • (54) 



The denominator in (54) takes into account the pole structure of these matrix elements 
whereas the function ^\™'."a n (/3i, • • • , b n ) carries the colour structure of this matrix ele- 
ment. The important point is that this function has neither poles or zeros in the physical 
sheet and satisfies a simplified system of monodromy equations which arise by substituting 
(54) into (46). 

Let us consider now the contribution of the four-particle intermediate states to the 
2-point functions. If the corresponding form factor was constant, this contribution would 
be of the form 



which corresponds in momentum space to a branch point at q = —16M Z of the form 
ln(g 2 + 16M 2 ). But since F m [ n vanishes at the threshold, neglecting other terms, we have 
an additional factor Yli<j Pfj m the integral. This leads to a behaviour ^—^ — , so that the 
branch cut gets softened to a function of the form (q 2 + 16M 2 ) 3 ln(q 2 + 16M 2 ). This effect 
actually gets stronger as we consider higher numbers of particles: the m-particle branch 
cut would be, on the grounds of phase space alone, of the form (q 2 + m 2 M 2 ) m /4-l ? 

but it 

actually gets softened to (q 2 + m 2 M 2 )™ 2 / A ~ l . 

As it is evident from the above discussion, this suppression of the higher particle states 
is not peculiar to the O(n) model but on the contrary we expect it to be a completely 
general feature for all interacting theories which have S(0) = — 1. From a pragmatic point 
of view, this observation is crucial to the success of the form factor bootstrap approach 
to the computation of correlation functions. In fact, although for integrable models the 
exact computation of the FF may be achieved with little effort on the basis of our previous 
considerations, one may wonder about their final usefulness. After all, to compute the 
correlation functions we have to sum over the infinite multiparticle FF, and these series 
are in general hard to analyse. Of course, for large values of Mr only the lowest term 



i<j 



cosh 




18 



will dominate the sum, the higher particle contributions being exponentially small. For 
small values of Mr however the correlators possess singularities and, on these scales, all 
numbers of particles contribute in principle to the sum. But if there is reason to believe 
that the higher particle terms are anyway suppressed, then the lowest contribution of the 
spectral representation becomes quite effective in approximating the correlation functions 
even in these ultraviolet regions. The validity of such approximation is of course linked to 
a sufficiently mild divergence of the correlation functions in the short distance scales. 

A first check of this phenomenon of suppression of higher FF contributions is given by 
looking at the second moment of the scaling function $i(p) 

U 2 = [ P 2 $i{ P )d 2 P . (55) 



The extra power in p 2 helps to smooth the ultraviolet behaviour of the correlation function 
and we may expect this integral to be almost saturated by the two-particle contribution. 
For the trace of the stress-energy tensor the two-particle form factor we take 1 



<0|6(0) \a,p 1 ;b,fc>= -2ttz <U M 2 F min (/?! - , (56) 

with -F m i n (/3) given by (53). Truncating the spectral representation of the scaling function 
<&l(p) to this term we have for the RHS of (55) 



d/3 

cosh 4 (3 ' 



^ I —r^\F{2P)\ 2 = 2.2241 , (57) 



to be compared with the exact result of U2, from (21), extracted by the c-theorem sum 
rule (in the limit n — > 0) 

20 

U 2 = — = 2.2222 . (58) 

Therefore we see that by keeping only the two-particle approximation of the correlator, we 
get an approximation to U 2 with a precision of about one part in 10 3 . However, the careful 
reader will notice that in this case the two-particle approximation results in an estimate 
of this quantity slightly larger than the exact value. Since the spectral representation of 
the scaling function <E>i(p) involves an infinite series of terms which are integrals over the 



We assume the mildest behaviour at infinity of this matrix element to fix it uniquely. 
Its normalization is easily determined by using the definition of the hamiltonian operator 

H = ^7r dx l T QQ (x° , x 1 ) and computing the matrix elements of both terms of this 
equation on (Pi,ai\ and |/32? a 2)- Notice however that our expression differs from that 
proposed in 20 which diverges at infinity much faster than ours. 
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modulus squared of its form factors, one may expect an underestimate of the exact value of 
U2 from the 2-particle approximation. Although this is what usually happens in ordinary 
QFT, the situation for the 0(n) (in the limit n — > 0) may be different. To see this, let us 
consider for instance the four-particle FF of the trace of the energy-momentum tensor 

(0\e\p 1 ,a;p 2 ,b;p 3 ,c;p 4 ;d) . (59) 

From group theory considerations this has the form 

A $ab$cd + B $adhc + C $achd > ( 60 ) 
where the invariant amplitudes A, £?, C (functions of Pi, . . . , /?4 and n) satisfy the functional 
and recursive equations of the form factors. For our considerations, their explicit form is 
not needed. In order to consider the four-particle contribution to the scaling function 
$l(p), we have to take the modulus squared of this form factor and sum on the colour 
indices. The final expression is given by 

(| A | 2 + I B | 2 + I C \ 2 )n 2 + (AB + BA + BC + CB + CA + AC)n , (61) 

which may be rewritten as 

(| A | 2 + I B | 2 + I C | 2 ) (n 2 -n) + | A + B + C | 2 n . (62) 

This is clearly positive for n > 1, and in this case we have the usual positive contribution 
of this form factor to the two-point function. But, for n — > 0, we need to pick up the linear 
term in n, i.e. the cross term in (61), which has no reason to be positive. As a matter of 
fact, for n = 1, which is the Ising model, we know it is negative, because in that case there 
is no 4-particle form factor so that A + B + C = 0. Hence, even though we have not pursued 
explicitly the calculation of the above amplitudes A, £?, C for the O(n) model, one should 
not be surprised to find that higher particle contributions will contribute negatively to the 
correlation functions. This effect becomes more evident if we consider the zero moment of 

*l(p) 

U = J $l(p)d 2 p . (63) 

This time there are no extra powers in p which suppress the ultraviolet behaviour of the 
correlation function, and higher form factors may expected to be important. As shown in 
the next section, Uq may be computed exactly through the TBA 

47T 2 

Uq = — = 13.1595 . (64) 
On the other hand, the two-particle approximation to this quantity is given by 

d 2) = 27T r°^ l F min(2t) | 2 = 14 3651 ^ (65) 

Jo coslr t 

which is close to the exact result, but slightly larger. 
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4. TBA AND THE RATIO C/D. 



4.1 The bulk free energy and the amplitude U. 



In this section, we discuss the calculation of the universal amplitude U in the extensive 
part of the free energy, defined by InZ ~ AnUM^, where M = is the mass of the 
0(n) model, in the limit n — > 0. This is performed using the thermodynamic Bethe ansatz 
(TBA). This method was first applied to integrable perturbed conformal field theories in a 
pioneering paper by Al. Zamolodchikov^, who also later^ applied it to the minimal models 
perturbed by the (1, 3) operator. Since our calculation is a rather simple generalization of 
his, we shall only sketch the argument, and refer the reader to Zamolodchikov's papers for 
the rather heavy details. Thus, in this section, we make no attempt to be self-contained. 

In the TBA approach, one considers the massive theory as a quantum field theory 
in 1 + 1 dimensions, at finite temperature R~^. Denoting the reduced free energy per 
unit length of this system by E(R), dimensional analysis implies that this has the form 
E(R) = (2n / R)F(r), where r = MR. As r — > 0, the scaling function has the form 
-F(O) + (9(r 2 ) + singular terms. The terms proportional to F(0) dominates in the high- 
temperature limit R — > and gives the central charge of the conformal field theory. The 
0(r 2 ) piece gives a term extensive in both the temperature and the spatial size of the 
system, and is therefore the extensive part of the free energy when we view the system as 
a classical system in two dimensions, which is what we wish to compute. 

The TBA proceeds by enumerating the allowed states in a large spatial box, labelled 
by the rapidities of the particles, consistent with the fact that the wave functions must 
change by a factor of the S'-matrix whenever two particles are exchanged. The fraction of 
these states which are actually occupied is then determined by minimizing the free energy. 
For a simple theory with only one type of excitation of mass M, the result is that 




(66) 



where 




(67) 



is the solution of the integral equation 



-MR cosh P + e(p) + ((f)* L) ((3) = 
21 



(68) 



where * denotes a convolution, and 4>{j3) = —i(d/d(3) In S((3). 

Zamolodchikov showed^ that, in the limit R — > 0, it is sufficient to consider the 'kink' 
solution of the simpler integral equation 

- e / 3 + e kink (/?) + (0*L kink )(/9) = , (69) 



and that 



I f C ^7 r kink / o\ 

F(r) ~ y L Vm \f3)e l3 df3 - ^ J KP) e^d(i + singular terms . (70) 

Furthermore, the integral in the second term may be evaluated easily by examining the 
behaviour of (69) as (3 — > — oo, and requiring that the 0(e@) terms cancel between the first 
and third terms. Thus if 4>([3) ~ Ae@ as (3 — > — oo, then the value of the integral in the 
second term of (70) is 2% /A. 

In Ref. 7, Zamolodchikov considered the case of the unitary minimal models, labelled 
by an integer k, which correspond to the critical continuum limit of the lattice RSOS 
models, perturbed by the operator (1,3). Although the O(n) model and the minimal 
model with n = 2 cos(7r/ (k + 1)) have the same value of the central charge, their operator 
content is certainly different. However, within the sector of operators of the type (1, 2s + 1), 
generated by repeated operator product expansions of (1,3) with itself, they are the same. 
Thus, to any order in perturbation theory, the correlation functions of the two theories 
should agree, and therefore so should such quantities as the mass gap and ground state 
energy. We assume that this will also hold non-perturbatively. In principle, these theories 
could differ at the non-perturbative level, but we see no physical reason for this to happen. 
Zamolodchikov argued that, although there is only one physical massive excitation in the 
RSOS model, in order to count correctly the states it is necessary to introduce pseudo- 
excitations which carry no energy but do enter into the TBA equations. For the kth 
minimal model these excitations correspond to the vertices of the Aj { _i Dynkin diagram, 
labelled by a = 1, . . . , k — 1. The physical particle corresponds to a = 1. The kink integral 
equation (69) becomes 

-5 al +e^ k ((3) + J2W*Lf- k )((3) = , (71) 



where is the incidence matrix of the diagram. Thus, by looking at the f3 — > — oo and us- 
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ing A = 2 we find the following system of equations for the integrals I a = j (clLjf nk / d(3)e ^d/3: 

I 2 = 7T 

h + h = o 

h + h = (72) 

Strictly speaking, these equations make sense only when k is an integer. However, we may 
formally find the solution for general k, which is 

.si n((j-fc)7r/ 2) 
sin(/c7r/2) 



h = * ■ ( 73 ) 



Since only a = 1 has mass, the term we want is proportional to I\ = tt cot(/c7r/2). Putting 
all the factors together, we find that 

F(r) = F(0) - — cot(/c7r/2) + power series in r 4 / ^ fc+2 ) . (74) 

87T 

Note that when k is an odd integer, the extensive part of the free energy vanishes, as 
found by Zamolodchikov: this is interpreted as being a consequence of (fractional) super- 
symmetry. For k = 2p the coefficient diverges, but, since the overall result is finite, there 
must be a cancellation of this leading term against one of the singular terms, resulting in 
logarithmic behaviour. This comes from 



k ^ 2p 8nk-2p\ ) k + 2\27iJ ' V ; 

in agreement with (2.22) of Ref. 7. 

However, we are interested in the limit n — > 0, corresponding to k — > 1. In that 
limit we then find for the free energy per unit area A~ l \nZ ~ (n/8)(k — 1)M 2 . Since 
dk/dn\ n= Q = 2/tt, we then have the final result that 

U = \ • (76) 
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4.2 The ratio C/D. 



We are now in a position to calculate the amplitude D defined in Sec. (2). From 
(21, 23, 76) we find that U 2 /U = 5/3tt 2 and so, from (24), 

D = A$ ■ (77) 

67T2 

Now consider the amplitude C, which, by (30), is proportional to the ratio of the second 
and zeroth moments of the scaling function $2- At the same level of approximation, 
we should neglect all m-particle intermediate states with m > 2. But, since the spin 
operator couples only to odd m, this implies that we need include only the single particle 
intermediate state, which is equivalent to assuming that 

.. , , f e iq -P d 2 q 

Note that this is not the same as assuming the Ornstein-Zernicke (free-particle) form for 
the correlation function, which would imply that r\ = and 7 = 2u. In writing (78) for 
the scaling function, we have already extracted in (30) the dependence coming from the 
exact value 7. From (78) we find the moment ratio V2/V0 ~ 4, so that finally, from (31), 
we obtain one of our main results 

C/D « 24 f F( f } « 13-70 . (79) 

5 r(§) 

It is difficult to estimate the errors in this calculation, which arise from the neglect of the 
3- and higher particle intermediate states in the spin-spin correlation function. We would 
expect the error to be larger for the zeroth moment Vq than for V2. On the grounds of 
phase space alone, one would expect the m-particle state to give a relative contribution 
of the order of J e~ mMr cfir ~ m~ 2 to the zeroth moment. This will be modified by a 
factor of m ~ TO ( m_ l)/4 ^ Q th e softening of the multi-particle branch points discussed in 
Sec. (3), leading to an expected error of a few percent in Vq. This estimate, however, does 
not take into account the further suppression of the small r region due to the small value 
of the magnetic scaling index 77. 
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4.3 Connection with the conformal limit. 



There is more information to be gained from the TBA. So far the energy operator 
£(r) has appeared only in the combination (2ir/v)(x c — x)£(r), as the trace 0(r) of the 
stress tensor. But if we fix on a definite normalization of £(r), for example, that given 
naturally by its conformal limit, then there is a definite relation between the temperature- 
like variable (x c — x) and the mass M. As shown by Al. Zamolodchikov 6 ' 7 , this may be 
determined by comparing the corrections to the TBA with the results of a perturbative 
analysis. This analysis has been carried through already 34 for the O(n) model for a slightly 
different situation than that considered above, namely, when the non-contractible loops 
which wind around the imaginary time direction carry a factor 2 rather than n. However, 
this change of boundary condition should not affect the relationship we wish to determine. 

With this boundary condition on the cylinder, the effective UV central charge is c = 1, 
rather than vanishing in the limit n — > 0. This may be traced to the existence of a negative 
dimension operator with scaling dimensions (Aq, Aq), such that 

c = c-12(A + A ) (80) 

so we see that Aq = Aq = — ^j, when n = 0. Noticing that these are the scaling dimensions 
of the operator <p2,3 m the Kac table, we may interpret this choice of boundary conditions as 
arising from the insertion of this operator at each end of the cylinder. Thus the perturbative 
expansion of the ground state energy begins as 35 

7r /27r\^ 3 
m) = -6R + R(Xc ' X) {-R) C f2,3),(2,3) + • • • > (81) 

where ^ ^ 3 ) is the structure constant of the three point function (£<P2,3 ( P2,3)- 
In our problem it is natural to normalize the energy operator so that as r — > 0, 

<5(r)5(0)> ~ n/r 4 / 3 , (82) 

so that, with periodic boundary conditions, the free energy will be order n as required. 
This means that we should take £ = y/nipi s, where ip\^ is normalized in the conventional 
fashion. As a result the structure constant is 

Cf , , , = lim v^d; (0 ^ = o — >- ; / = 0.302248 . (83) 

(2,3), (2,3) k _ >1 V V ) (2,3), (2,3) 16yr 3 T(4/3) V ' 

Note that this has a finite limit as n — > 0. 
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Defining the scaling function c(MR) = —(6R/n)E(R) we therefore have 

c(MR) = 1 - 12(2tt)" 1 / 3 (0.302248) (x c - x) # 4 / 3 + • • • . (84) 



On the other hand, in Ref. 34 c(MR) has been determined as the numerical solution of 
the TBA equations, and the first terms of its expansion are given by 

c(MR) = 1 - 0.4454536 (MR) 4 / 3 + ■■■ . (85) 

Comparing these, we obtain the number which links the physical mass M to the scale of 
temperature, given our choice of normalization of £ 

(x c -x) = kM A /' s , (86) 

where k 0.226630. Moreover, using our parametrization M = ^q 1 ^ — x/x c ) 3 / 4 , we also 
see that 

x c = kZ 4/3 . (87) 

This relation will be used in section 5 to fix the asymptotic behaviour of the structure 
factor. 



5. THE STRUCTURE FACTOR OF SELF-AVOIDING LOOPS 



An important measure of self-avoiding walks, which is, in principle, accessible by light 
scattering experiments, is the structure factor of iV-step loops 

2 



S N (q) = N- 



iq-n 



Min-rj) 



) ■ 



(88) 



Sj^(q) is a positive definite quantity that plays the role of generating function for all 
moments of the energy-energy correlator. In the continuum limit we have 



n 



J2N PN S N (q)x N - 2 = a [ d 2 xe^ x (£(x)£(0)) , 
N J 



(89) 



and Sjy(q) may be expressed as 



S N (q) ^ NF(y) = N 1 - *- + . . . 



(90) 
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where y is the scaling variable given in terms of the mean square radius of gyration of 
loops 

y 2 = q 2 (R 2 g ) N . (91) 

As we are going to show, the scaling function F(y) may be analyzed in great detail by 
means of two different methods, the form factor approach and the conformal field theory. 
By the former we are able to arrive at quite an accurate estimate of the scaling function 
F(y) (at least for moderate values of the scaling variable y), whereas, by the latter we may 
extract its asymptotic behaviour. 



5.1 Structure Factor in the Two-Particle Approximation 



A precise determination of F(y) may be obtained in terms of the two-particle form 
factor as follows. Expand initially both terms in (89) in power series of q 



1 



oo . 

nJ2 N PNS N ,2 P Q 2p x N - 2 = a J2^ J d 2 x(iq ■ x) 2p (S(x)S(O)} . 



(92) 



N,p p=Q 



(the odd terms in the series vanish by parity). The key point is that we may compute 
exactly the first two coefficients of the series by exploiting the TBA and the c-theorem sum 
rule. As the other terms are related to the higher moments of the energy-energy correlation 
function 

J d 2 x\x\ 2p {S(x)S(0)} , (93) 

they will be approximated with quite high accuracy by retaining only the two-particle 
contribution. The reason is that the higher powers \x\ 2p emphasize more the large dis- 
tance scales of the correlation function, and on the contrary suppress its short distance 
singularity. Hence they improve substantially the ultraviolet convergence of the integral 
and the spectral representation of the correlator (S(x)S(0)) is effectively truncated to the 
two-particle term 1 . 

Let us then write the right hand side of (92) as 

a (h ~ \h + h) , (94) 



As a matter of fact, we have already seen in sect. 3 that the two-particle contribution 
provides an accurate estimate for the second moment. Hence, we expect the precision 
obtained by the two-particle contribution to increase for the higher moments. 
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where 

h = J d 2 x(£(x)£(0)) 

I 2 = J d 2 x(q-x) 2 (£(x)£(0)} (g5) 

oo 

As discussed in Sec. (2), I\ and I2 are exactly obtained in terms of the extensive part of 
the free energy and by the c-theorem: 

2v(2v- 1) o rr 

^3 is related to the higher moments of (£(x)£(0)), and, computed in the two-particle 
approximation, it is given by 

h= 2n(x c -x)^ 2 ( - 1) m^ l2p ' (97) 

where 

r^^¥^L ■ (98) 



/ (cosht) 2 + 2 P 

The numerical evaluation of the first is reported in Table I. They vanish asymptotically 



as 



l2 " P -» 4 r( P+ l) _( ' ' ( ' 

where S(0) is a constant defined in the appendix. 

To obtain S]^(q) from the RHS of (92) we need to express Mas fo^l-x/xc) 3 / 4 and 
to take the inverse Laplace transform. The result is 

$N,2 5 T(N) 



£d 



s Nfi [ } 2^ 2*p r (at - 1) r (3£+il) ^° 
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3(p-l)\ (100) 



The non-universal factor £o as in (100) may be eliminated by using both the scaling form 
of Sw(q) ((90)) and the corresponding definition of (Rg)]\f obtained by 5/v,2/£/V,0 

in agreement with (77). With this substitution, the scaling function F(y) is given by 

(2 00 \ 
i-y + EWH > ( 102 ) 

where 

Their values for p < 20 are reported in Table I. The resulting series is convergent for all 
complex y, and in particular yields a highly precise determination of the scaling function 
F(y) for small values of y. Although the higher order coefficients become more and more 
precisely determined by the 2-particle approximation, this does not mean that this ap- 
proximation captures the correct large y behaviour of the scaling function. Indeed, it is 
straightforward to see that, within the 2-particle approximation, this is determined by the 
large (3 behaviour of F m \ n ((3), which leads to a power law A more accurate way of 

finding the large y behaviour is through conformal field theory. 



5.2 Asymptotic Behaviour of the Structure Factor from CFT 



In order to analyze the asymptotic behaviour of F(y) in the limit y — > 00, let us first 
write the RHS of (92) as 

MQ) = r dttUt)(S{t/q)Sm . (104) 

<T JO 

(Jg(t) is the Bessel function). In the limit q — > 00 only the short distance scales of 
the correlation function are important and therefore (£{t/q)£(0)) is determined by the 
operator product expansion: 

(£(r)£(0)) ~ r~ 4 / 3 + C (£(0)) r~ 2 / 3 + . . . . (105) 
r— >0 

(£(0)) is the vacuum expectation value of the energy operator determined by TBA 
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whereas C is the structure constant of the CFT algebra. To compute it, we have to use the 
formulas of Ref. 36 and take into account the fact that we have chosen to normalize the 
2-point function of the energy operator S(r) to n rather than unity. Therefore C is given 
by the limit 

' V (2 



C = lim v^) C^ M (k) = 6 ( j . (107) 



Inserting (105) into (104), the large q behaviour of A(q) is given by 

A(q) ~ A^ + A^ + ... 

q^oo 



2tt a J\ 3 a M 2 CJ 2 1 



(108) 



+ 



where 



^2/3 32n(x c -x) g 4/3 

roc r 
Ji = / dtr 1 ' 2, J (t) = 2 1 / 3 — 
Jo r 

j 2 = r d tt^j ( t ) = 

io r (I) 



(109) 



To obtain the scaling function F(y) we still have to perform an inverse Laplace transform 
on A(q). However the leading order A^ 2 / 3 ) of A(q) does not contain any dependence on 
(x c — x) and therefore in the asymptotic expansion of F(y), 

the coefficient a vanishes identically! Hence the scaling function F{y) decreases at infinity 
faster than expected by a power-counting argument. As we will discuss at the end of this 
section, the validity of this result is not restricted to two-dimensional self-avoiding loops 
but, on the contrary, is quite general. In two dimensions, however, we may use the exact 
value of the structure constant C to extract the universal coefficient b in (110). In fact, 
expressing M as Co" 1 (l — x/x c ) 3 / 4 and making an inverse Laplace transform on A^l^ we 
have 

N PN S N W, ~ _^C« -\ c ^f 9 -V» . (ill) 
Normalizing this quantity to the corresponding expression for q — > 0, i.e. 
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(in agreement with (7,76),) and using the relation (R 2 ) N ~ 5£qN 3 / 2 /6ir 3 / 2 from (77), 

4/3 

we see that the coefficient b in (110) depends on the combination k = x c ^q , which was 
computed in Sec. (4.3). The appearance of this quantity in this context is not surprising, 
since we are considering the ultraviolet behaviour of S(q), which is normalized at q = 0. 
Putting all this together, we obtain the final result 

F[V) ^ + ' ' ' ■ (U3) 
where b is the universal constant 



(114) 

9.65065 x 10~ 3 . 




(2n)l 



Returning to the question of the leading coefficient a of F{y) which vanishes, it is 
easy to show that this result is quite general and may be simply justified by a physical 
argument. In fact, in d dimensions the leading behaviour of the generating function of 
Sjq{q) decreases as q d ^l v ? but with a corresponding coefficient which is independent of 
(x c — x). The sub-leading term goes in any dimension as q^l v but the coefficient in front 
depends on the contrary on (x c — x). Therefore, the scaling function F(y) always decreases 
at infinity as y^l v . This is what one expects for an object whose fractal dimension is \ jv. 
It is easy to see that the same power law also rules the large q behaviour of the structure 
factor of a linear polymer whose generating function is given by 

e i? ( r - r ')(s(0)^)^')s(ri)> c AA'Ai . (115) 



Since this is a connected correlation function, in the limit q — > oo the leading term comes 
from the £ term in the OPE of £(r)£(r') rather than from the unit operator. Hence the 
decrease of F(y) as y~ l l v is consistent with the idea that linear and loop polymers have 
the same fractal dimension on small scales. 



6. AREA OF LOOPS AND THE 
CURRENT-CURRENT CORRELATION FUNCTION. 



In this section we show how the amplitude governing the mean square area of self- 
avoiding closed loops is related to the second moment of a current-current correlation 
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function in the 0(n) model. A summary of this calculation was given in Ref. 10. Consider 
a loop of N steps. It has been argued that the mean area (a) ^ of such loops should 
behave asymptotically as EN 2 ", and, more generally, that the moments (a p )]\[ of the area 
should behave as E^N 2pv , where the amplitude combinations E^/D p are expected to 
be universal. In this paper we focus on the case p = 2. 

Any given loop may be assigned two possible orientations. Consider a given orientation, 
and let jft^(r) be a vector of unit magnitude on the link r in the direction of the orientation. 
Then the signed area of a given loop is 

a = ^V^W • ( 116 ) 

r 

Of course, this quantity averages to zero when summed over orientations, but its mean 
square is non-zero. For a given loop 

« 2 = 1 E (^*V)*V) - r,r' v J l f{r)J^{r')) . (117) 

r,r' 

Using the results ^ r «^i at ( r ) = and ^ r r v J^{r) = 0, which follow from the fact that J v 
is conserved, (117) may be rewritten as 

« 2 = JE(-^- r0 2 jj*(r) #V) + (r M - r^){r v - r' v )J^{r)J^{r')) . (118) 

r,r' 

This is for one loop. If we now average over all possible oriented loops, weighted by a 
factor x N for each loop, we find that 

2nN s ^p N (a 2 ) N x N 
N 

= \ E ( " U r - r') 2 (jl at (r)J^(r')) + (r, - ^)(r v - r' v ){J^\r)J^\r'))) 
r,r' 

(119) 

where the correlation functions on the right hand side are evaluated in a complex O(n) 
lattice model, with hamiltonian = —x J2 nn s*(r)s a (r'), in the limit n — > 0. The current 
is the one which generates the U(l) transformations s a e za s a . An explicit expression 
for J^ at (r) was given in this model in Ref. 10. 

(119) is now ready for taking the continuum limit. (Doing this earlier on (117) leads 
to erroneous results.) To do this, we simply let ^ r J^ at (r) — > f J^{r)d 2 r, so that the right 
hand side of (119) becomes 

\A J{-\r 2 5^ + r^r v ){J^r)J v {Q))d 2 r . (120) 
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Since is conserved, by rotational symmetry the correlation function on the right hand 
side has the form 

(J M (r)J„(0)> = {d^dv - 5^d 2 )G(r) , (121) 

where G is a scalar. Substituting this into (120) and integrating by parts, the right hand 
side is simply A J Gd 2 r. On the other hand, using the same method, one may show that 
/ r 2 (J M (r)J M (0))rf 2 r = -4/ Gd 2 r. Thus 

2nJ2PN(a 2 ) N x N = ~la Q I r 2 (J^{r)J^))d 2 r . (122) 
N J 

Now, since the current has unit scaling dimension, its two-point function has the scaling 
form 

(J M (r) J M (0)> = -nrHj(r/i) . (123) 

Moreover, the normalization of this current is fixed by the requirement that the ends of 
a self-avoiding walk, which in the complex O(n) model are represented by insertions of 
s\ and si, are respectively unit sources and sinks for J^. Therefore the scaling function 
<E»j, like that of the correlation function of the trace of the stress tensor, is completely 
universal with no metric factors. It follows that the right hand side of (122) has the form 
lna Uj${l - x/x c y 2 \ where Uj = J p 2 §j{p)d 2 p. Thus 

PN(a 2 ) N ~\aaQt 2 Q U JT ^-p N . (124) 

Using p N ~ BN~ 2v - l ii N and the relation (20) for BD, we then find that 

Finally, using the exact relation (77) for D, we obtain 

(a 2 ) N 48tt 3 



(&)% ~ 25 



Uj . (126) 



Now let us see how we might compute the current- current correlation function, and 
hence Uj, through the form factor approach. Since the current is an O(n) singlet, it couples 
only to the even particle sector. The simplest non-trivial form factor is therefore 

(0\J^x 1 m,a 1 ,+;P 2 ,a 2 ,-) = e m ^ h ^+^ h . 

(127) 
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In the above ± denotes the U(l) charge of the particle, and aj its O(n) colour index. It is 
simpler to work with the light-cone combinations = Jq±J\. Then, by O(n) symmetry, 
Lorentz covariance and current conservation, 

(01^(0)1/51, ai, +; p 2 , a 2 , -) = <W 2 (e ±A - e ± ^)F(/9 12 ) . (128) 

It is straightforward to check that this satisfies the conservation condition q + J~ + q~ J + = 
0, with = M(e ± ^ 1 +e ziz ^ 2 ). In fact the Watson equations for this form factor then show 
that the simplest solution is to take F([3) = dF m [ n (P), where d is a normalization constant, 
and -F m j n was defined in Sec. (3). The normalization is fixed by crossing (128) to find 

(Pl,ai,-\J ± (x 1 )\p2,a2,-) = -rfVa 2 (e ±A +e ± ^)F mi Jz7r-/3 12 )e" jM ( sinh / ? - sinh ^)- 1 . 

(129) 

Thus the matrix element of the total U(l) charge is 



(/?l,a b -| Jj°(x 1 )dx 1 \p2,a 2 ,-) 



= -27T d 5 ai a 2 (cosh Pi + cosh/?2)-F m in(«7r - /?i2)<!> (M(sinh Pi - sinh/^)) . 

(130) 

The left hand side is simply — 27r5 ai a 2 ^(Pl ~ Therefore, using the fact that 

F m j n (z7r) = 1, we find 

d=\M . (131) 
In the 2-particle approximation, then, the current- current correlation function is 

(WIPDl = ^ / |i^(eA-eft)(e-A-e-A)|F mk (ft 2 )|^-MKco S h A+ c„ 8 h & ) 

(132) 

The factor of n comes from the sum over colour indices in the intermediate state. Note 
that there is no factor of 1/2! since the particles have opposite U(l) charge and so are not 
identical. 

This leads to the estimate for Uj in this approximation 

Uj3 = -^i ^ dl> ■ (133) 

The integral may be evaluated numerically using the forms for F m [ n discussed in Sec. (3), to 
give U Jj2 = 0.61506701. The resulting estimate for the ratio E^/D 2 is 36.62. This figure 
is unreasonably high. Estimates of the mean area 32 lead to (E^/D) 2 « 6.39, and we 
would not expect this to be significantly less than E^/D 2 . Moreover, since the maximal 
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area of tiR is reached for a circular loop, we would not expect this amplitude ratio to 
exceed 7r 2 . In fact, recently Guttmann 37 has estimated a value for for the square 
lattice which gives pa 6.72 a much more reasonable value. 

With hindsight it is not hard to find the reason for the failure of the 2-particle ap- 
proximation in this case. If one examines the region of integration which is contributing 
the dominant part to the integral in (133) one finds that it comes from \(3\ ~ 2 — 5, which 
corresponds to a center of mass energy between 7.5M and 150M, which is far above the 
4-particle threshold at AM. Thus, under these circumstances, there is no justification 
for ignoring the 4-particle (and perhaps higher) intermediate states. The fact that the 
2-particle contribution overestimates the result is not surprising in the light of our result 
for Uq, where we showed that higher intermediate states can give negative contributions 
in the n — > limit. There appear to be two reasons for this failure. First, the contribution 
near the 2-particle threshold is suppressed by a factor (1 — cosh 2/3), whose origin may 
be traced to current conservation. Second, the {J^{r)J v {^))) correlation function behaves 
like r in the ultraviolet limit, and is therefore much more singular than the spin-spin 
or energy-energy correlation functions. Thus it is much harder to approximate this be- 
haviour keeping only the two-particle state, and the ultraviolet region will give a much 
more important contribution to the second moment. 



7. RADIUS OF GYRATION OF SELF-AVOIDING WALKS. 



In the previous sections, we have shown how various amplitudes of self-avoiding walks 
and loops are related to moments of two-point correlation functions in the O(n) model, 
which may be estimated using the form factor approach. In this section we consider another 
such quantity, which is, however, related to a higher correlation function. The evaluation 
of this therefore provides a further test of the form factor method. 

Consider and A-step self-avoiding walk from the origin to the point r. If labels a 
site visited by the walk, then the squared radius of gyration is given by 

2N 2 Rl = J2(r l -r J ) 2 . (134) 

hi 



Following the same line of argument as in Sec. (2), if cjy the total number of A-step walks, 
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the generating function cj^N^{R^) j^x N is proportional to the moment 

{n -r 2 ) 2 (s(0)£(ri)S(r 2 )s(r))d 2 rid 2 r 2 d 2 r . (135) 



/< 



The calculation of this four-point function using the form factors would be very cumber- 
some: depending on the time-ordering of the points, different kinds of intermediate states 
would arise, and, in addition, it would be necessary to know at least the three-particle 
form factors of the spin operator s, which we have avoided so far in this paper. 

However, the radius of gyration is in fact related to other measures by the Cardy-Saleur 
formula 38 (as corrected by Caracciolo et a/. 39 ), which is a consequence of the c-theorem: 

^(Rl) N -2(R 2 m ) N + ±(R 2 e ) N = , (136) 

where (Rm)N is the mean square distance of a monomer (a site visited by the walk) from 
one end of the walk, and, as in Sec. (2), (R 2 )n is the mean square end-to-end distance. 
Thus, if we can compute one universal ratio of these measures, we may find the other. In 
fact, (Rm) n is easier, because it involves calculating only a three-point amplitude. 

In fact, let us consider the integral 

J = ! (r 1 + r 2 ) 2 {£{Q)s{r l )s(r 2 ))d 2 r l d 2 r 2 . (137) 



It will become clear below why we choose this particular moment. Since {r\ + r 2 )^ = 
2r\ + 27*2 — (r\ — r 2 )^, this is related to generating functions in the following way: 

J = 4 c N N(R 2 m ) N x N - 1 - c N N{R 2 e ) N x N ~ l . (138) 
N N 

Let us now see how to evaluate J using the form factor approach. In general we have to 
sum over all possible time orderings of the points 0, r\ and r 2 , with a different expression 
in terms of the form factors in each region. In addition, the regions where t\ and t 2 are 
on the same side of the origin will involve, even in the simplest approximation, knowledge 
of the form factors (/3i/3 2 \s\/3^). We may avoid this by the following trick. Let v\ = R + p 
and r 2 = R — p. If we imagine doing the p integral first, the result will depend only on 
\R\, so we may choose the time component of R to be zero. This means that the time 
ordering will always be sSs and not Ess or ss£. Let us therefore write p = (t, x) and 
assume that the intermediate states are saturated by their lowest contribution, which, in 
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this time ordering, is only the 1-particle state. Then 
J-^J 4R 2 (£(0)s(R- p)s(R + p))d 2 Rd 2 p 

= 64tt fVdfl r dx Hdt f^f ^(olslfoe-imR-^h-MtcoBhfk 
JO J-oo Jo J 2?r J 2n 

■ (A | £ I (3 2 > e lM {R+x) sinh p2 ~ Mt cosh 02 (fc | s 1 0) 

(139) 

In writing this we have suppressed the colour indices, which are all set equal to a fixed 
value by the insertion of the spin operator, rather than being summed over. The integrals 
over x and t are straightforward (it is for this reason that we choose the particular moment 
defined by J.) The form factors of the energy operator £ are simply proportional to those 
of the trace of the stress tensor, 6, given by (56), so that 

(p\£\ - (3) = i M 2 u(x c - x)- l F min {m + 2(3) . (140) 

This gives 

roo r jo — 2iMRsmh(3 

J = 327ns(x c -x)- 1 \(0\s\l)\ 2 R 3 dR -f iF min {m + 2(3) . (141) 

Jo J 2n cosh z /? 



Before performing the integral over R it is advantageous to deform the contour so that 
f3 — > (3 — in/2, with the new (3 contour lying just above the real axis. This is possible 
because of the analyticity properties of F m j n . The R integral is then straightforward, and 
we find 

J tt§v{xc-x)~ l M- A \(id\s\l)\ 2 I , (142) 

where 

Jc cosh 4 (3 sinh 2 (3 

The contour runs just above the pole at (3 = 0. The contribution from this small semicircle 
gives 7rS(0), where the number H(0) is evaluated exactly in the Appendix. The remaining 
principal value integral may be evaluated numerically using the form for F m [ n given there. 
The result is that I w ttS(O) - 0.456 « 3.432. 

Now, in the same approximation, 



J2cn(RI)nx N = [ r 2 (s(r)s(0)}d 2 r 
N J 



J l<°l-l 1 >| 23, /(M^ I = Mf - 4 |<°l-l 1 ) 



(144) 
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so, for consistency with (27), we should take (up to a non-universal normalization which 
will cancel) 

|(0|s|l)| 2 = M r i . (145) 
Thus, in the one-particle approximation, 

Y J C N {R 2 e ) N x N = ^ + \l-x/x c y 1 - 2v , (146) 
N 

so that 

cnV&n ~ ^ff^rf N ■ ("7) 

Now, using the same normalization, the desired combination of generating functions 
(138) is 

4j2c N N(R 2 m ) N x N -J2c N N(R 2 e ) N x N = - x/xc)-^' 1 , (148) 

N N 
so that 

AT7+2i/+l 

c N N{A(R 2 m ) N - (R 2 e ) N ) ~ 6i/Jg T{l + 2v + 1 f N > ( 149 ) 

and finally 

ARla)N 1 = 6^J r( 7 + 2z/) 6^7 
<^1>AT 8 r( 7 + 2z,+ l) 8( 7 + 21/) ' 1 J 

Substituting our estimate for 7, we then find (R^) n / (R 2 ) n ~ 0.420, which corresponds 
to I \R 2 )n{R 2 ) N ~ 0.126. The best numerical estimates for these ratios, from Monte Carlo 
simulations, are 0.43962 ± 0.00033 and 0.14029 ± 0.00012 39 . 

The approximately 5% error in the approximation of the first ratio is what might be 
expected from the earlier results of this paper. The particular moment chosen in (137) 
does not damp the short-distance behaviour of the correlation function as r\ or become 
small. On the other hand, the singularities in this correlation function are less severe than 
those in the energy-energy correlation function, and in that case, an estimate of Uq led to 
an error of less than 10%. Unfortunately, the error in our estimate gets magnified when 
we use the Cardy-Saleur formula (136) to estimate the other ratio. 

Clearly it is possible to improve this calculation at the expense of considerably more 
analytic work. However, our result does show that the utility of the form factor approach 
is not restricted to the computation of two-point functions. 
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8. DISCUSSION. 



In this paper we have attempted to show that the study of two-dimensional field theory 
applied to critical behaviour has advanced to a point where it is now capable of yielding 
not only exact values for critical exponents but also other universal amplitude ratios and 
scaling functions which describe the approach to the critical point. We have chosen the 
self-avoiding walk problem as a non-trivial example where precise numerical comparisons 
may be made with the results of series expansions. 

Our main results are the estimate for the ratio C/D, given in (79); the calculation 
of the scaling function of the structure factor for closed loops, whose the coefficients of 
whose power series expansion are tabulated in Table I, and whose asymptotic behaviour is 
given in (113, 114); and the amplitude ratio for the radius of gyration of open self-avoiding 
walks, given in Sec. (7). 

The S'-matrix and form factor approach to correlation functions, while not exact in the 
sense that the sum over intermediate states must be truncated, turns out to be remarkably 
useful in practice. In most cases, only a few-particle intermediate states need be retained 
to give a precision of one part in 10 3 . This is related to the softening of the branch cuts 
by the interactions. This effect is peculiar to two dimensions and is a consequence of 
two-particle unitarity. Consider, for example, a scalar field theory with a g<p^ interaction. 
Two-particle phase space is the imaginary part of the simple bubble diagram 

d 2 k 1 



/ 



(k 2 + M 2 )((q - k) 2 + M 2 ) MvV+ 4M 2 



(151) 



However, in the interacting theory this is modified by the iteration of the bubble diagram 
so that it becomes 

1/g (152) 



MvV + 4M 2 1 + g/(M^/q 2 + 4M 2 

as q 2 — > — 4M 2 . Thus, no matter how small g, the two-particle branch cut is always 
softened. A similar effect is responsible for the behaviour S — > — 1 of the two-body S- matrix 
at threshold. The Watson equations also show that these two results are manifestations 
of the same phenomenon. For higher particle states, the action of this effect in each 
two-particle channel separately further suppresses their contribution. 

For the self-avoiding walk problem, this effect has a simple physical interpretation: the 
repulsive interaction between different parts of the walk suppresses those configurations 
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where the walk repeatedly winds back on itself. It is these which correspond to the multi- 
particle intermediate states in the form factor approach. 

While this method works well in most of the examples considered in this paper, it 
fails in others. The degree of success appears to be directly related to the softness of the 
short-distance singularities in the correlation function, which is given by the conformal 
dimensions of the operators involved. For example, for the spin-spin correlation function, 
which has a singularity of the form r~5/ 24 , even the simple one-particle approximation does 
very well, even for the zeroth moment. For the energy-energy correlation function, which 
behaves like r -4 /^, the 2-particle approximation works to a level of 10% for the zeroth 
moment, and to 1 part in 10^ for the second. The TBA and the c-theorem give these 
moments exactly, in any case. In the three-point function (sSs), which has a singularity 
r -37/48^ \ eve \ Q f accuracy is about 5%. However, for the current- current correlation 
function, which behaves like r , the two-particle approximation fails completely for the 
low moments. It would therefore be useful to have another method of interpolation between 
the infrared and ultraviolet behaviours of such correlation functions. 

The possibility of computing higher-point correlation functions also suggests that these 
methods may give other interesting geometric information on the sizes and shapes of self- 
avoiding walks and loops, for example, their asphericity. 
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APPENDIX 

In this appendix we discuss some useful formulas for the F m i n (f3) in the limit n — > 0. 
For large values of (3, F m [ n (f3) behaves as 

Fmin(P) ~ . (153) 
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Its analytic structure may be read off from its infinite product representation 



F min (P) = sinh^£(/?) 



(154) 



where 



so?) = n 



r(z + i + |f)r(z + 2 + ^ 



2^ 



r z i + 



r 2 (/ + | + I 



r(/ + i)r(/ + 2) 



2(1+1) 



(155) 



and /3 = iir — (3. An interesting quantity is the value of H(/?) at the origin 



H(0) = exp (^C(3)) = 1-23756 



(156) 



(where £(3) is the Riemann function). 

For numerical calculations, a useful formula is given by a mixed representation of 

2 



N-l 

m = n 



x exp 



1 - ' m 



1 + 



mi 
i+i 



i + 



/9/27T 

z+2 



2 



(157) 



"I 



^ dxsmh^fq(N,x)e~ 2N7TX 
a; sinh 2 7rx cosh ^ 



sin 



/9a; 



where 



q(N,x) = 



(N + 1)(N + 2) _ ^ + 2) e _ 27r , + (jVCjV + 1) e - 47 rx 



(158) 



2 v 2 

In fact, the rate of convergence of the integral may be improved substantially by increasing 
the value of N. 
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Table I 



V 


%2p 


B 2p 





2.28627 


1 


1 


3.53980 x 10" 1 


-0.5 


2 


1.65158 x lO" 1 


1.17363 x lO" 1 


3 


1.00385 x 10" 1 


-1.65011 x 10~ 2 


4 


6.92166 x 10~ 2 


1.58450 x 10~ 3 


5 


5.14233 x 10~ 2 


-1.12325 x 10~ 4 


6 


4.01424 x 10~ 2 


6.18847 x 10~ 6 


7 


3.24617 x 10~ 2 


-2.74819 x 10~ 7 


8 


2.69541 x 10~ 2 


1.01085 x 10~ 8 


9 


2.28462 x 10~ 2 


-3.14532 x 10~ 10 


10 


1.96857 x 10~ 2 


8.41962 x 10~ 12 


11 


1.71925 x 10~ 2 


-1.96581 x 10~ 13 


12 


1.51847 x 10~ 2 


4.04936 x 10~ 15 


13 


1.35396 x 10~ 2 


-7.43064 x 10~ 17 


14 


1.21715 x 10~ 2 


1.22477 x 10~ 18 


15 


1.10193 x 10~ 2 


-1.82639 x 10~ 20 


16 


1.00382 x 10~ 2 


2.47954 x 10" 22 


17 


9.19444 x 10~ 3 


-3.08162 x 10~ 24 


18 


8.46266 x 10~ 3 


3.52340 x 10~ 26 


19 


7.82306 x 10" 3 


-3.72244 x 10" 28 


20 


7.26018 x 10~ 3 


3.64834 x 10~ 30 



Table I. The values of the first twenty integrals X^p defined in (98), and the coefficients B 2p of the structure 
factor scaling function (102), expanded in powers of its argument y — qR. 
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